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We propose a new form of inhomogeneous phases consisting of out-of-phase staggered flux domains 
separated by diagonal charged domain walls centered either on bonds or on sites. Remarkably, such 
domain flux phases are spin-rotationally symmetric and exhibit cone-like quasiparticle dispersion 
near the Fermi energy, as well as incommensurate order of orbital currents. Such features are 
consistent with the pseudogap behavior and the diagonal stripes observed experimentally in lightly 
doped cuprates. A renormalized mean field theory shows that these solutions with coexisting charge 
modulation and charge currents are competitive ground state candidates within the t—J model. 

PACS numbers: 74.72.-h, 71.45.Lr, 74.20.Mn, 75.40.Mg 



I. INTRODUCTION 

Among numerous new ideas and concepts that have 
been put forward to explain the unusual properties the 
high temperature superconductors (HTS), which go be- 
yond the conventional Fermi liquid theory^ the staggered 
flux (SF) phase 2 attracts much attention as a candi- 
date for the pseudogap normal phase of the underdoped 
cuprates. 3 Such a state is characterized by a checkerboard 
pattern of plaquette currents circulating clockwise and 
anticlockwise on two different sublattices so that the cor- 
responding flux flowing through each plaquette alternates 
in sign. 

On the one hand, using the SU(2) gauge invariance 
of the Heisenberg model one can show that at half-filling 
the SF phase is equivalent to the d-wave superconducting 
wave function^ which has correctly reproduced several 
key experimental properties of the HTS. 5 Moreover, its 
Gutzwiller-projected energy is in a very good agreement 
with the best estimate for the ground-state energy of the 
two-dimensional undoped Heisenberg antifcrromagnet. 1 
On the other hand, even though a finite doping removes 
this degeneracy and stabilizes d-wave superconductivity 
in the ground stated the SF phase is the lowest-energy 
Gutzwiller-projected nonsuperconducting state that has 
been constructed so far^ and its energy spectrum re- 
mains similar to the d-wave superconductor. Signa- 
tures of the SF pattern in the current-current correla- 
tion have been seen in the Gutzwiller-projected d-wave 
superconducting phased and in the exact ground-state 
wave-function of the t—J models It has also been pro- 
posed that the hidden d-density wave (DDW) order of 



the doped SF phase could be the origin of the mysterious 
pseudogap behavior J£ Finally, it has been shown that 
under some circumstances the SF phase can coexist with 
d-wave superconductivity in the underdoped regime^ 

However, the physics of the hole-doped cuprates seems 
to be even more involved as the competition between 
the superexchange interaction which stabilizes the anti- 
ferromagnetic (AF) long-range order in the parent Mott 
insulator, and the kinetic energy of doped holes, might 
lead to the formation of stripe phases with hole-rich re- 
gions and locally suppressed magnetic order, which was 
suggested in early Hartree-Fock studies. 1 - In a stripe 
phase two neighboring AF domains are separated by a 
one-dimensional domain wall (DW), where a phase shift 
of 7T occurs in the AF order parameter. Later on, ex- 
perimental confirmation of the stripe phases has trig- 
gered a large number of studies devoted to their prop- 
erties within a number of methods which go beyond the 
Hartree-Fock approach*^ Moreover, even though static 
charge and spin orders have only been observed in layered 
cuprates, e.g., in Lai.6-a;Ndo.4Sr a Cu04 (Nd-LSCO) (see 
Ref. [Tib and La2- a; Ba 2 ;Cu04 (see Ref. fl5h . while in bi- 
layered YBa2Cu306+<5 (YBCO) only a stripe-like charge 
order and incommensurate spin fluctuations have been 
reported^ stripe phases quickly joined the list of can- 
didates for the pseudogap phase in the cuprates as they 
are compatible with many experimental results^ 

Although numerical simulations of microscopic models 
of correlated fermions, such as the t—J model (see later), 
are especially difficult, various signatures consistent with 
(i) DDW states, and (ii) stripe phases have been de- 
tected. In particular, the emergence of strong staggered 
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current correlations under doping the Mott insulator has 
been reported in exact diagonalizations by Leung^ and 
attributed to the formation of spin bipolarons. 18 These 
findings are consistent with an early observation of stag- 
gered spin chiralitjsi^ since charge degrees of freedom 
strongly couple to spin scalar chirality. Interestingly, 
spin chirality /charge currents seem to compete with hole 
pairing, 20 and this issue requires a further careful consid- 
eration. Simultaneously with those findings, the obser- 
vation of stripes and checkerboard patterns (which also 
include some form of charge ordering) has also been con- 
firmed by density matrix renormalization group (DMRG) 
computations for some boundary conditions i^i 

We also note that an exotic SF phase with long-range 
orbital current order at half-filling (in contrast to the 
fully projected SF phase, see Ref. y|) was stabilized in var- 
ious extended Hubbard-like models (which include some 
form of charge fluctuations not present in the simpler 
model discussed above) within ladder— or bilayer— ge- 
ometries. It was also shown that such a long-range DDW 
order could survive with the emergence of stripe-like fea- 
tures under doping^ 

Unfortunately, even though stripe phases seem to play 
important role in the physics of HTS, it is still not clear 
how the stripes are connected, as a competing state, to 
d-wave superconductivity. Therefore, in this paper we 
introduce a new class of wave functions with compos- 
ite order in a form of filled domain flux (FDF) phases, 
with one doped hole per one DW atom. In addition 
to capturing essential properties of the SF phases, the 
FDF structure accounts for the incommensurate diag- 
onal spin peaks observed in lightly (x < 0.06) doped 
La 2 _ a; Sr a; Cu04 (LSCO)^ and Nd-LSCO^ Thus, our 
phase should allow one to obtain a smooth transition 
from the insulating state at half-filling to the d-wave su- 
perconductor above a critical doping x c , with a concomi- 
tant change of the DW orientation into vertical stripes 
just at x c , as observed experimentally in LSCOi^l The 
existence of such phases is suggested by recent variational 
Monte-Carlo calculations which show an instability of the 
SF states towards phase separation^ and we argue that 
self-organization into flux domains separated by DWs is 
generic in the doped t— J model. Most pronounced fea- 
tures of these phases shown in Fig. [T](a,b) are: (i) doped 
holes self-organize into diagonal DWs, (ii) DWs separate 
weakly doped SF domains with a smoothly modulated 
magnitude of the flux within them, (iii) DWs introduce 
a phase shift of 7r in the flux phase and the SF domains 
alternate, and finally (iv) in contrast to the so-called com- 
mensurate flux (CF) phases, the total flux vanishes, and 
therefore no asymmetry of the magnetic response is ex- 
pected when reversing the direction of an applied mag- 
netic field. In fact, these FDF phases have strong similar- 
ities with the solution obtained in Ref. [28| using uniform 
(i.e., site independent) Gutzwiller factors. 

The paper is organized as follows. The t-J model and 
its treatment in the Gutzwiller approximation are intro- 
duced in Sec. |TTJ The properties of locally stable domain 



flux phases with either bond-centered or site-centered do- 
main walls are presented in Sec. IIII1 The paper is con- 
cluded in Sec. HY]by pointing out certain possibilities of 
experimental verification of the suggested type of order 
and by a short summary of main results. 

II. MODEL AND FORMALISM 

We consider the t-J model^ 

« = - E *« + h.c.) + ^ • S,-, (1) 

(ij),o (ij) 

which is believed to describe the physics of the HTS. 5 
Here the summations include each bond (ij) only once. 
Next, the local constraints that restrict the hopping pro- 
cesses oc cl^Cjp to the subspace with no doubly occu- 
pied sites are replaced by statistical Gutzwiller weights^ - 
while decoupling in the particle-hole channel yields the 
following mean field (MF) Hamiltonian, 

Hmf = - E *«^i( c Ls> + h - c -) - V E Uia 

{ij),cr ia 

~ 1 J E 9ij(Xji4a c ^ + h - c - ~ \Xij\ 2 ), (2) 

with the self-consistency conditions for the bond-order 
parameters 

Xji = (Cj a c ia ). (3) 

In principle, simultaneous decoupling in the particle- 
particle channel is also possible^ but since we are inter- 
ested in the diagonal DWs similar to the ones observed 
in the underdoped LSCO famil y 25 ! 26 we focus here on 
nonsuperconducting solutions. In particular we choose 
x = 1/16, one of the magic doping fractions at which 
low-tempcraturc in-plane resistivity of LSCO is weakly 
enhanced suggesting a tendency towards charge order^ 
Here, to allow for small non-uniform charge modulations, 
the Gutzwiller weights have been expressed in terms of 
local doped hole densities 

»M = l-5>Ua> (4) 
cr 

as follows: 33 

9ij = y/^Zj, 9ij = (2 - Zi)(2 - Zj ), (5) 

with z% = 2nhi/(l + nhi)- For simplicity, results shown 
below correspond to nearest neighbor hopping tij = t 
only. 34 Thanks to developing an efficient reciprocal space 
scheme by making use of the symmetry— the calcula- 
tions were carried out on a large 256 x 256 cluster at low 
temperature f3J — 500, which eliminates the finite size 
effects. 
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FIG. 1: (color online) Spatial modulation of the hole density nhi (circles), bond-order parameter \ij (lines with arrows 
indicating the direction of charge currents), and flux $□ defined by Eq. (|6} (positive/negative flux indicated by symbols +/— ) 
distribution found in two FDF phases at hole doping x = 1/8 and t/J = 3. Circle diameters are proportional to the doped 
hole densities; widths of bond lines connecting them are proportional to the magnitudes of the bond-order parameters Xiji 
while the magnitude of flux flowing through each plaquette is represented by the size of +/— symbol. Two distinct phases are: 
(a) bond- centered FDF phase with a vanishing current (dashed lines) at the DW bonds; (b) site-centered FDF phase with a 
vanishing flux (indicated by 0) at the DW plaquettes. Panel (c) shows the self-consistent CF phase (t — 0) characterized by 
the uniform fictitious flux $□ = ^(1 — x), as well as by homogeneous charge distribution. 



Our starting point is the CF phase, a wave function 
which, away from half-filling, displays remarkable com- 
mensurability effects at special fillings and fulfills the self- 
consistency condition at t = Indeed, in the limit 
xt/J — > 0, the magnetic (superexchange) energy in the 
CF phase exhibits a minimum when the fictitious flux 
(in unit of the flux quantum), flowing through each pla- 
quette and defined by a sum over the four bonds of the 
plaquette 

= E %> ( 6 ) 

(ij)eO 

where 0y is the phase of Xij > follows exactly the filling 
fraction, i.e., — i(l — x). In this case, Hamiltonian 
(J2j) reduces to the Hofstadter Hamiltonian describing the 
motion of an electron in a uniform magnetic flux assumed 
to be rational <!>□ = p/q^ Therefore, the peculiar prop- 
erty of the superexchange energy follows from the CF 
phase band structure with q bands and the Fermi level 
lying in the largest gap above the pth subband. As a 
result, the modulus of the bond-order parameter x-ij ©, 
the spin correlation and the hole density are all spatially 
uniform [see Fig. HJc)]. However, infinitesimally small 
xt/J selects a special arrangement of the phases {0^} 
so as to optimize the kinetic energy term oc cos ®*j 
and should produce an inhomogeneous structure.— 

Within this class of singlet (nonmagnetic) wave func- 
tions, competing with possible inhomogeneous solutions 
(see later) , the uniform SF phase also offers a very good 



compromise between the magnetic (Ej) and kinetic (E t ) 
energy. For small t and x, the kinetic energy is mini- 
mized (within the MF approach) when all phases of \ij 
are set to a constant 0y = ±7r/4, corresponding to al- 
ternating fluxes $□ = ±0.5 (SF phase). Increasing xt/J 
gradually reduces |$rj| and drives the system towards a 
Fermi liquid state (with real Xij) m a continuous way. 

III. DOMAIN FLUX PHASES 

Starting with initial parameters corresponding to a 
uniform CF phase, the self-consistent procedure leads 
to new FDF solutions which could explain a diagonal 
spin modulation observed experimentally in the insu- 
lating regime of LSCO^ and Nd-LSCO^ usually in- 
terpreted in terms of diagonal stripes, even though no 
signatures of any charge modulation were observed yet. 
This conjecture is also supported by the recent neutron 
scattering studies of the Ni impurity effect on the diago- 
nal incommensurability in LSCOj^ Indeed, doping by Ni 
quickly suppresses the incommensurability and restores 
the Neel state. This indicates a strong effect on hole lo- 
calization and thus favors the presence of charge stripes 
with mobile holes rather than the spiral order with local- 
ized hole spins. 

Interestingly, we found two types of topologically dif- 
ferent but nearly degenerate solutions which both have 
the same size of the unit cell (see Fig. [T|): (i) a bond- 
centered FDF phase, very similar to the original CF one, 
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where each DW is characterized by a zero current stair- 
case and by a maximum of the hole density spread over 
the related bonds [Fig. []Ja)], as well as (ii) a site- centered 
FDF phase, where the DWs are characterized by zero flux 
plaquettes ordered along a diagonal line and by a maxi- 
mum of the hole density centered at two of their corner 
sites [Fig. [TJ b)] . Apart from local doped hole densities 
{nhi}, bond quantities are needed for a full characteriza- 
tion of both phases (here we use a short-hand notation) : 
— the spin correlation 



&i ~ 29'li+x\Xi,i+x\ 2 , 

the bond charge hopping 

the charge current 

U = ^9i,i+x Im {Xi,i+x}, 
as well as the modulated flux 



2,2+2; : 



(7) 



(8) 



(9) 



(10) 



with a phase factor (—1 compensating the modula- 

tion of the flux within a single domain of the SF phase. 
Typical profiles of the above defined observables at low 
doping are depicted in Fig. [2j 

The stability of the FDF phases originates from a sub- 
tle competition between the magnetic Ej and kinetic en- 
ergies E t . Let us first focus on the t/J—>0 limit where 
the site-centered SF phase is stable and very competitive 
(among the nonmagnetic states), in contrast to the bond- 
centered one. This extreme case corresponds to the lo- 
calization of doped holes at DWs and the superexchange 
energy in the SF domains is best optimized. Indeed, by 
expelling holes from the SF domains one reinforces locally 
the AF correlations with a concomitant reduction of both 
bond charge and current correlations. On the contrary, 
due to a large hole density, both these tendencies are re- 
versed around the DWs. However, increasing t/J leads 
to a much broader charge spatial distribution in the unit 
cell as a larger fraction of holes enters the SF domains 
(see Fig. [2]). Nevertheless, both FDF phases remain com- 
petitive even in the regime of large (realistic) values of 
t/J ~ 3 due to: (i) enhanced short-range AF correla- 
tions deep in the SF domains (Si ~ —0.33 compared to 
S ~ —0.28 in the uniform phase), where the fictitious 
flux approaches the special value $ = h (local minimum 
of the kinetic energy in the limit xt/J — > 0), and (ii) 
strongly enhanced bond charge accumulated around the 
DWs, typically three times larger than that in the SF 
phase, due to both amplification of the gh factors and 
reduced (vanishing) fictitious flux flowing through the 
bond-centered (site-centered) plaquettes at the DWs. 

Of particular interest is whether one can also stabilize 
within the present formalism the so-called half-filled do- 
main flux (HDF) phases, analogous to half-filled stripes 
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FIG. 2: (color online) (a,e) Hole density nhi (0|, (b,f) spin 
correlation Si 0, (c,g) bond charge Ti ©, and (d,h) modu- 
lated flux c&Tri (|10[) in the bond- centered (left) and site-centered 
(right) FDF phases at x = 1/16 for: t/J — 1 (triangles), and 
t/J = 3 (squares). For comparison, circles depict the related 
t/J — > solutions: the CF phase with uniform fictitious flux 
$ = 15/32 (left) and a two-domain |$| = | SF phase (right). 



with one hole per two atoms in a DW as observed in the 
cuprates around x = l/8i 14 ' 15 On the one hand, both 
self-consistent bond- and site- centered HDF phases found 
at x = 1/16 and t/J = 3 have a somewhat higher total 
energy per site [F ~ —1.03 J) than those obtained for 
both degenerate FDF ones (F ~ —1.07 J), and for the 
uniform SF phase (—1.09 J). However, Table Q] shows 
that all domain flux phases become very competitive at 
x = 1/8, not only with respect to the SF phase but also 
with respect to a recently proposed nonuniform 4x4 
superstructure^ Note also that while the FDF phases 
optimize mainly Ej, the HDF ones are characterized by 
rather low E t . Therefore, we predict that large t/J rather 
favors the domain flux phases with partially filled DWs. 
We argue that quantum fluctuations are likely to stabi- 
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TABLE I: Kinetic energy per hole Eh (in units of i), and 
kinetic energy Et, magnetic energy Ej, free energy F (all per 
site in units of J) for the locally stable phases: bond- centered 
HDF(l) site-centered HDF(2), 4x4 checkerboard, FDF, and 
SF one, as found at hole doping x — 1/8 and t/J = 3. FDF(l) 
and FDF(2) phases are fully degenerate. The lowest energy 
increments are given in bold characters. 



phase 


E h 


E t 


Ej 


F 


HDF(l) 


-2.7856 


1.0446 


-0.4028 


-1.4474 


HDF(2) 


-2.7843 


-1.0441 


-0.4026 


-1.4467 


4x4 


-2.7128 


-1.0173 


-0.4348 


-1.4521 


FDF 


-2.7067 


-1.0150 


-0.4418 


-1.4568 


SF 


-2.7587 


-1.0345 


-0.4246 


-1.4591 



lize them, in analogy to the half- filled stripe phases^ or 
to the fully projected 4x4 checkerboard wave function 
which was recently shown to be more stable than the 
uniform SF phased This suggests that other inhomoge- 
neous solutions might be stable as well. Unfortunately, a 
direct comparison of our singlet wave functions to the 
original (magnetic) stripe phases^ 2 - is not possible yet 
since both are described within two entirely different for- 
malisms. Hence further studies using more sophisticated 
methods (like projected wave functions as in Ref. l38h are 
needed. 

An experimental support of the FDF phases fol- 
lows from angle-resolved photoemission (ARPES) exper- 
iments on lightly doped LSCO that show a strongly sup- 
pressed spectral weight near the pseudogapped X — 
(ir, 0) and Y = (0, n) points, and a quasiparticle band 
crossing the Fermi energy /i along the nodal T — M direc- 
tion, with M = (7r,7r).— Both features are qualitatively 
reproduced in the FDF phases - the electronic bands are 
almost dispersionless along the X—Y direction, and a gap 
opens at w = /j (Fig. [3]) , indicating that transport across 
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FIG. 3: (color online) Electronic structure of the site-centered 
FDF phase (solid lines) and SF phase (dashed lines) along 
the main directions of the Brillouin zone for x = 1/16 and 
t/J — 3. Inset shows a pseudogap between the FDF bands 
along the X — Y direction near the Fermi energy \i (thin 
dashed line). 



the DWs is suppressed. However, the most salient feature 
of the electronic structure in FDF phases is a relativis- 
tic cone-like dispersion around the S = (tt/2, ir/2) point. 
Indeed, massless Dirac excitations are at the heart of the 
quantum electrodynamics in (2+1) dimensions (QED3) 
theory of pseudogap in the cuprates.— This feature is 
also found in the SF phase, but for the uniform flux and 
hole distribution it occurs away from the Fermi energy fi. 
The shape of the electronic structure in the FDF phase 
depends on the actual value of t/J. Firstly, a strong lo- 
calization of holes at DWs in the limit t/ J — > pushes 
the top of the lower band cone well below /1. Secondly, 
finite t weakens the stripe order so that the gap between 
the lower and upper band at the S point is reduced. A 
further increase of t pushes some lower band states above 
\i enabling transport along the DWs. 



IV. DISCUSSION AND SUMMARY 

For possible experimental verification of the present 
proposal it is important to realize that orbital cur- 
rents of the domain flux phase give rise to weak mag- 
netic fields (that should be experimentally distinguish- 
able from the copper spins). Muon spin rotation (/iSR) 
technique is an extremely sensitive local probe especially 
suited to study small modulations of local fields. Earlier 
estimations^ give 10 to 100 Gauss corresponding roughly 
to 0.03 to 0.25 lib in cuprates. In fact, incommensu- 
rate order in the LSCO family seen in neutron scattering 
measurements ; 25 ' 26 (with a large but finite correlation 
length) might be attributed, at least partly, to the ex- 
istence of orbital moments. Finally, note that although 
the phases considered here do not break SU(2) symme- 
try and do not exhibit AF long range order, on general 
principle they can still sustain AF correlations on large 
distances (i.e., beyond nearest neighbor sites) between 
copper spins. 

In summary, we have introduced and investigated a 
new class of flux phases that unify the remarkable prop- 
erties of the SF uniform phase with the incommensu- 
rate magnetic correlations established in the underdoped 
cuprates. Bond- and site-centered FDF phases are nearly 
degenerate which indicates strong fluctuations which are 
expected to be amplified, either for increasing t/J or for 
increasing doping x. As these phases are only marginally 
unstable at the MF level, they might be stabilized by 
quantum effects and explain the low temperature physics 
of the cuprates in the low doping regime, where a pseu- 
dogap phase forms at higher temperature. Therefore, 
the solutions presented here could be viewed as a low- 
temperature instability of the nearby DDW pseudogap 
phase (stable at higher temperature but below T*) in the 
same way as the "ordinary" stripe phases could be seen 
as an instability of the nearby doped AF Neel state at 
infinitesimal x. Therefore, our proposal calls for a search 
of experimental signatures of domain flux phases in the 
underdoped cuprates, especially in the LSCO family. 
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